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Abstract 
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included in the Higgs mass calculation. A link to MicrOMEGAs allows to compute 
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up to the GUT scale. 
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1 Introduction 



The Next to Minimal Supersymmetric Standard Model (NMSSM) [1-9] provides a 
very elegant solution to the \x problem of the MSSM via the introduction of a singlet 
superfield S. For the simplest possible scale invariant form of the superpotential, the 
scalar component of S acquires naturally a vacuum expectation value of the order of the 
SUSY breaking scale, giving rise to a value of /i of order the electroweak scale. Hence the 
NMSSM is the simplest supersymmetric extension of the standard model in which the 
fundamental Lagrangian contains just SUSY breaking terms but no other parameters of 
the order of the electroweak scale. 

In addition, the NMSSM renders the "little fine tuning problem" of the MSSM, orig- 
inating from the non-observation of a neutral CP-even Higgs boson at LEP II, less se- 
vere [2]. 

As in the MSSM, the phenomenology of the NMSSM depends on a certain number 
of parameters (mostly soft SUSY breaking parameters) that cannot be predicted from an 
underlying theory at present. It is then useful to have computer codes that compute 
physically relevant quantities as Higgs and sparticle masses, couplings, decay widths etc. 
as functions of the initial parameters in the Lagrangian. Such codes allow to investigate 
which regions in parameter space are in conflict with present constraints on physics beyond 
the standard model and, most importantly, which regions in parameter space can be tested 
in future experiments and/or astrophysical measurements. 

In the MSSM, corresponding computer codes are HDECAY [10], FeynHiggs [11], 
Isajet [12], SoftSusy [13], MicrOMEGAs [14], Suspect [15], Spheno [16], CPSUPERH [17], 
SDECAY [18] and DARKSUSY [19]. In the NMSSM, the only available code at present 
is NMHDECAY [20]. 

In the present paper we describe the improvements performed in the version 2.0 of 
NMHDECAY. (In the meantime version 2.1 is available, whose features are described in 
the file README on the web page [20]. In the cases where they differ from the ones of 
version 2.0 they are described below.) First we recall the features of the previous version 
of NMHDECAY, version 1.1. Starting from a set of (low energy) parameters it performs 
the following tasks: 

• It computes the masses and couplings of all physical states in the Higgs sector and 
in the chargino and neutralino sectors. 
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• It computes the branching ratios into two particle final states (quarks and leptons, 
all possible combinations of gauge and Higgs bosons, charginos, neutralinos, but 
not decays into squarks and sleptons) of all 6 Higgs particles of the NMSSM. Three 
body decays via WW* and ZZ* are computed as in HDECAY [10], but no four 
body decays are taken into account. 

• It checks whether the Higgs masses and couplings violate any bounds from negative 
Higgs searches at LEP, including many quite unconventional channels that are rele- 
vant for the NMSSM Higgs sector. It also checks the bound on the invisible Z width 
(possibly violated for light neutralinos). In addition, NMHDECAY 1.1 checks the 
LEP bounds on the lightest chargino and on neutralino pair production. 

• It checks whether the running Yukawa couplings A, k, h t or encounter a Landau 
singularity below the GUT scale. 

• Finally, NMHDECAY 1.1 checks whether the physical minimum (with all vevs non- 
zero) of the scalar potential is deeper than the local unphysical minima with van- 
ishing (H u ), (H d ) or (S). 

The improvements in the versions 2.0+ are as follows: 

1. Further radiative corrections are added in the Higgs sector, in order to improve 
the precision of the Higgs masses calculations. In addition, all squark and slepton 
masses (and mixing angles for the third generation) are computed. 

2. Branching ratios of all Higgs states into squarks and sleptons are computed, and 
squark and slepton loops are included in the Higgs decays to two gluons and two 
photons. 

3. Experimental constraints from LEP and Tevatron on squark, gluino and slepton 
masses are checked. 

4. The dark matter relic density can be computed, via a call of a NMSSM version of 
MicrOMEGAs (that is provided on the same web site). 

5. The branching ratio BR(b — > 57) is computed to lowest order. 
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6. In the version 2.1, the RGEs for the soft /susy/ breaking terms can be integrated 
up to Mqut- 

7. A Makefile for optimal compilation is provided. 

The conventions concerning the Lagrangian of the model are the same as in version 1.1: 
The superpotential W is given by 

W — h t Q ■ H u f c R - h b Q ■ H d B R - h T L ■ H d L c R + XSH U ■ H d + ^kS 3 . (1.1) 

(Hereafter, hatted capital letters denote superfields, and unhatted capital letters the cor- 
responding (complex) scalar components.) The SU(2) doublets are 




Products of two SU(2) doublets are defined as, e.g., 

H u -H d = H+H d - H° u H° d . (1.3) 

For the soft SUSY breaking terms we take 

-£ S oft = m 2 Hu \H u \ 2 + m 2 Hd \H d \ 2 + ml\S\ 2 + m 2 Q \Q 2 \ + m 2 T \T 2 \ 
+m 2 B \B 2 R \ + m\\L 2 \ + m 2 \L 2 R \ 
+(h t A t Q ■ H U T R - h b A b Q ■ H d B c R - h T A T L ■ H d L c R 

+XA X H u -H d S+l kA r S 3 + h.c.) . (1.4) 

The resulting mass matrices and couplings can be found in the appendix of ref. [20]. The 
conventions are also listed within the FORTRAN code as comments at the beginning of 
each corresponding subroutine. 

The input parameters relevant for the Higgs sector of the NMSSM (at tree level) are 

A, k, A\, A K , tan/5 = (H u ) / (H d ) , fi cS = X (S) . (1.5) 

As in the case of version 1.1, it is possible to a) use input and output formats according 
to the SUSY Les Houches Accord (SLHA) conventions [21] (with, however, a modified 
switch (65 instead of 23) for /z e fr), b) use a privately defined input and output format 
(scan), that allows to scan over a user defined range of the input parameters (1.5). 
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In the next section, we describe the improvements of versions 2.0 and 2.1. Apart 
from the additional radiative corrections, we discuss the precise meaning (renormalization 
scale) of the input parameters. In section 3, we describe in detail how the different versions 
NMHDECAY can be installed, compiled, and linked to MicrOMEGAs. We conclude with 
a short outlook. 

2 Improvements in NMHDECAY 2.0 and 2.1 

2.1 Radiative corrections to the Higgs masses 

Concerning radiative corrections induced by (s)top and (s)bottom loops, we adopt the 
same philosophy as in version 1.1: First we compute the running Yukawa couplings ht, 
hb and A as well as the Higgs vevs (and tan/3) at the scale of the stop and sbottom 
masses, taking into account effects of order h^ b times potentially large logarithms. (The 
input parameters A (and hence /i e fr), k, A x , A k are defined at the scale M§, which we can 
identify with mf op within the present approximation.) Then we add the one loop radiative 
corrections ~ h^ b to the effective potential, taking the full dependence on stop/sbottom 
masses and mixing angles into account. The advantage of performing this calculation at 
the scale of the stop/sbottom masses is that the remaining dominant two loop corrections 
~ /iL ~ h-t/b as are relatively simple. These latter corrections include now also the 
dependence on hb (in contrast to the version 1.1) in order to cover the large tan/5 regime. 
Also, as in version 1.1, due to this procedure several non leading two loop effects (related 
to squark mass splittings and mixings) are automatically included, as in ref. [22] for the 
MSSM. 

New in the version 2.0 are the corrections of the order ~ 9 2 h^ b to tne CP even 
Higgs boson masses (where g denotes the electroweak gauge couplings), induced by the 
stop/sbottom D term couplings, whose dependence on the stop/sbottom masses is com- 
puted beyond the leading logarithmic approximation. 

Concerning the logarithmic one loop corrections of the order ~ g A (to the mass of 
the lighter doublet like CP even state), we distinguish now the masses of the different 
squarks/sleptons of the different generations (assuming the first two generations to be 
degenerate). New in the version 2.0 are the logarithmic one loop corrections to fourth 
order in the NMSSM specific Yukawa couplings A and k, that have only recently been 
computed in ref. [23]. 
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Finally, in version 2.0 the corrections of the order ~ 9 2 h 2 / b to the Higgs pole masses 
(from contributions ~ h 2 j b to the Higgs self energies as in ref. [24]) are included. 

The dominant sources of uncertainty on the mass of the lighter doublet like neu- 
tral Higgs boson mass are thus non logarithmic contributions of the orders g 4 , g 2 \ 2 , A 4 , 
g 2 n 2 , k a and \ 2 k 2 , and two loop contributions beyond the dominant double logarithms 

~ h t/b> h t/b as - 

2.2 Higgs decays 

In addition to the Higgs decays considered in the version 1.1, the version 2.0 of NMHDE- 
CAY includes now the possible two body decays of all CP even, CP odd and charged Higgs 
bosons into all squarks and sleptons. Also, squark and slepton loop contributions to the 
radiatively induced decays of the CP even and CP odd Higgs bosons into two photons 
and two gluons are included. 

The corresponding code is essentially as in HDECAY [10], but with the more compli- 
cated Higgs self couplings and mixings in the NMSSM. The leading logarithmic corrections 
from top and bottom quark loops to the Higgs self couplings are included, and the Higgs 
squark couplings are scaled up to a scale Q 2 corresponding to the squark masses which 
takes care of large logarithmic radiative corrections ~ h 2 ^ b \n(Mg quark /M^), neglecting 
terms ~ h 2 /b ln(M 2 H igg J M 2 quark ). 

2.3 Sparticle masses 

The masses and mixing angles of the two charginos and five (in the NMSSM) neutralinos 
were already calculated in the version 1.1 of NMHDECAY in the subroutines CHARGINO 
and NEUTRALINO. In the actual version 2.1, one loop radiative corrections to the neu- 
tralino and chargino mass matrices are included as in section 4.2 in [24]. 

In the subroutine MSFERM in the version 2.0+, also the slepton masses are computed, 
and the squark mixing angles and pole masses are calculated including the one loop a s 
corrections^ . 

The subroutine GLUINO in the versions 2.0+ includes a computation of the gluino 
pole mass to the order a s . 

In the actual version 2.1 we assume that the input values of the soft SUSY breaking 
§We thank S. Kraml for contributions to the corresponding codes 
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terms are given at a SUSY scale Q 2 = (2Mq + Mfj + Mf,)/4, where M Q , M v and M D are 
the running squark masses of the first two generations. If desired, this scale can also be 
set by the user. 

Thus, in the versions 2.0+ of NMHDECAY the complete sparticle spectrum is com- 
puted. Already in the version 1.1, the masses of the two charginos and the masses and 
couplings to the Z boson of the five neutralinos were compared to LEP constraints from 
direct searches and constraints on the invisible Z width. Now, in addition, NMHDE- 
CAY tests the squark and gluino masses against constraints from the Tevatron [25] and 
LEP [26]. (The Tevatron constraints are those used by the LEPSUSY Working group [26].) 
As usual, NMHDECAY issues a warning in case where any of the present constraints is 
violated. NMHDECAY is thus quite unique in testing the complete Higgs and sparticle 
spectrum against constraints from accelerator experiments. 

2.4 Dark matter relic density 

The dark matter relic density in the NMSSM has recently been studied in ref. [27] (for 
previous investigations, see refs. [28]). For this purpose, NMHDECAY 1.1 was used to 
compute the Higgs and sparticle spectrum, which was then fed into a new version of 
MicrOMEGAs extended to the NMSSM. MicrOMEGAs calculates all the relevant cross- 
sections for the lightest neutralino annihilation and coannihilation. It then solves the 
density evolution equation numerically, without using the freeze-out approximation, and 
computes the relic density of the lightest neutralino. 

In the version 2.0 of NMHDECAY, the dark matter relic density can be computed for 
any choice of input parameters, by setting a simple flag in the input file, through a link 
to this NMSSM version of MicrOMEGAs. The details on how this link has to be installed 
will be given in the next section. 

In case the corresponding flag is on, the computed amount of dark matter is compared 
to constraints from WMAP [29] (.0945 < Qh 2 < .1287), and a warning is issued in 
case the result is too large or too small. (Clearly, a warning corresponding to a too small 
LSP relic density can be ignored, if one is ready to assume additional contributions to the 
dark matter.) 
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2.5 BR(6 -> s 7 ) 

In the version 2.0 of NMHDECAY the branching ratio BR(6 — > S7) is computed to lowest 
order in the subroutine BSC Contributions from charged Higgs and chargino/squark 
loops are included, which are the same as in the MSSM. The theoretical error is about 
6% for tan f3 < 3, but increases rapidly with increasing tan (3. Hence the result is only a 
rough estimate for tan (5 > 5. 

2.6 The RGEs 

A new feature of version 2.1 is the (possible) evolution of all NMSSM parameters up to 
the GUT scale. Three different subroutines are present: 

• the subroutine RGES integrates the RGEs for the gauge and Yukawa couplings up 
to the GUT scale (defined by the matching of the U(1) Y and £77(2) gauge couplings 
g 1 and g 2 ). It is always called and serves to compute M GUT as well as the Yukawa 
couplings at the GUT scale in order to check the absence of a Landau singularity 
(by requiring any Yukawa coupling to be less than Att). Two loop /3 functions are 
used, and the conversion from the MS to the DR scheme as well as all possible 
threshold effects between Mz and the SUSY scale are taken into account. 

• the subroutine RGESOFT integrates the RGEs for the couplings and the soft terms 
up to the GUT scale (that is now known from a previous call of RGES). This 
subroutine is always called in the version SLHA, but in the version SCAN it is called 
only if the output format "long" is chosen. (It is not useful, but time consuming, 
for large scans.) At present, terms ~ A 2 , k 2 in the two loop coefficients of the j3 
functions for the soft terms are missing (this will be improved in the near future). 

• the subroutine RGESOFTINV integrates the RGEs for the couplings and the soft 
terms from the GUT scale down to the SUSY scale Q. It is called only if RGESOFT is 
called. Unless modified, it uses the outputs of the subroutines RGES and RGESOFT 
for the values of all parameters at the GUT scale and leads thus to no new results for 
the parameters at the SUSY scale. Unless modified, it produces no output and acts 
just as a dummy. However, the user can easily modify the values of the (or some) 
parameters at the GUT scale (by choosing INGUT=1, after which the subroutine 
has to be re-compiled), and generate an output for the resulting parameters at the 
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SUSY scale. These can subsequently be used as input, and allow to generate sets 
of parameters that have desirable properties (as partial universality) at the GUT 
scale. Clearly, the Higgs masses squared at the SUSY scale cannot be put in directly, 
and can be modified only indirectly by varying other parameters like tan/3, /j, e s, A, 
k, A x or A K . 

3 How to install NMHDECAY 

Two versions of NMHDECAY are available: 

1. nmhdecay_slha uses an input file and produces output files that are suitable gener- 
alizations of the SLHA conventions [21]. It is configured for studying the properties 
of one user-defined point in parameter space. 

2. nmhdecay_scan employs privately defined input and output files. It allows to scan 
over parts of or all of the NMSSM parameters A, k, // e ff — \{S), tan/9, A\ and A K . 

Both programs are based on one single Fortran code (nmhdecay_slha.f or nmhdecay_scan.f) 
contained in the compressed directory N MH DECAY. tgz that can be downloaded from the 
web page http://www.th.u-psud.fr/NMHDECAY/nmhdecay.html. 

The directory NMHDECAY contains also a Makefile as well as input files for both 
versions and test output files. Once the compressed tar file NMHDECAY. tgz is downloaded, 
one should type: 

tar -zxvf NMHDECAY. tgz 
cd NMHDECAY 
. / make 

The command make will create 2 executable files, nmhdecay_slha and nmhdecay_scan. If 
one wishes to compile just one of them, it suffices to type ./make nmhdecay_slha or 
./make nmhdecay_scan. Both codes need data files in order to check Higgs, sfermion 
and gluino mass bounds from LEP and Tevatron. These data files are contained in the 
compressed tar file EXPCON.tgz that can be downloaded from the same web page. To 
uncompress it, simply type 

tar -zxvf EXPCON.tgz 



9 



which generates a directory EXPCON. By default in NMHDECAY, the path to the files 
containing the constraints is ../EXPCON so that the directory EXPCON should be located 
at the same place as the directory NMHDECAY. However, if one wishes to place it some- 
where else, it is possible to define a corresponding environment variable EXPCON_PATH: 

a) From a C-Shell, type 

setenv EXPCON_PATH f ull_path_to_EXPCON 

For a permanent setting add the same line to the file .cshrc. 

b) From a Bash Shell, type 

export EXPCON_PATH=full_path_to_EXPCON 

For a permanent setting add the same line to the file .bashrc. 

Finally, for those wishing to compute the dark matter relic density in the NMSSM, a 
version of MicrOMEGAs_1.3 extended to the NMSSM [27] is contained in the compressed 
tar file micromegas_l.3_nmssm.tgz, also available on our web page. In order to uncompress 
and to compile the C code omg.c and the MicrOMEGAs library, one should type: 

tar -zxvf micromegas_l.3_nmssm.tgz 
cd micromegas_l . 3_nmssm 
. /micro_make omg.c 

This will generate the executable file omg that can be called from NMHDECAY. It is pos- 
sible to put the directory micromegas_1.3_nmssm wherever one wants, provided one defines 
the environment variable MICROMG_PATH as above. Otherwise micromegas_1.3_nmssm 
has to be put in the same directory that contains NMHDECAY and EXPCON (ie the 
default value for MICROMG_PATH is ../micromegas_1.3_nmssm). However, if one wants 
to move the directory micromegas_1.3_nmssm (and consequently change the variable Ml- 
CROMG_PATH) after the MicrOMEGAs library has been compiled, one needs to recompile 
it. In order to do this, one has first to remove the generated files using the script clean 
(included in the directory micromegas_1.3_nmssm). The sequence of commands is then: 

mv micromegas_l . 3_nmssm new_path 
setenv MICROMG.PATH new_path (C.Shell) 
or 
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export MICROMG_PATH=new_path (Bash-Shell) 
cd new_path 
. /clean 

. /micro_make omg.c 

We now outline the particular features of the two versions of NMHDECAY. 

3.1 NMHDECAY_SLHA 

The program nmhdecay_slha uses the input file slhainp.dat, aversion of which (slhainp.dat. test) 
is contained in the directory NMHDECAY. This sample file appears in Table 1 below. Sev- 
eral comments on its contents are in order: 

a) BLOCK MODSEL contains the switch 3 (corresponding to the choice of the model) 
with value 1, as attributed to the NMSSM in ref. [21]. Any other choice would cause 
the program to stop. We have also added the switch 9 corresponding to a flag for the 
call to MicrOMEGAs: if the flag is set to 0, the relic density is not calculated (and Mi- 
crOMEGAs does not need to be compiled). If the flag is set to 1, nmhdecay_slha calls the 
executable file omg which returns the relic density of the lightest neutralino if it is the LSP. 

b) BLOCK SMINPUTS contains important Standard Model parameters. 

1. The first is the inverse electromagnetic coupling constant a~^, now at the scale 
Q = M z (in contrast to the version 1.1 where it was taken at the scale Q — 0). 

2. Second, since Higgs vevs and couplings are defined in terms of Mz, My/ and Gf, an 
on shell scheme is used implicitly, and the values of Mz, My/ and Gp are required. 
Mz, Gp (and a s (M z )) are defined in BLOCK SMINPUTS, whereas the numerical 
value of My/ is defined in the subroutine INITIALIZE. This subroutine assigns de- 
fault values to all parameters and reads the files with the experimental constraints 
(contained in the directory EXPCON). 

3. Third, as part of this block the running b quark mass m^m^), the top quark pole 
mass and m T are read in. 

4. In addition, nmhdecay_slha needs the strange quark mass m s {MS) and the charmed 
quark (pole) mass (taken as 190 MeV (at Q = 1 GeV) and 1.42 GeV, respectively, 
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as in HDECAY [10]), as well as the CKM matrix elements V us , V c b and V u b- The 
numerical values of these five parameters are defined in the subroutine INITIALIZE. 
(For convenience, they are printed out in the output file spectr.dat, see below.) 

c) BLOCK MINPAR contains the optional switch with the input value for the renormal- 
ization scale Q at which the squark and slepton masses, trilinear couplings and gaugino 
masses are defined. For Q = (or no switch 0) this renormalization scale is computed in- 
ternally from the average of the first generation squark masses. The switch 3 corresponds 
to the input value for tan/?. 

d) BLOCK EXTPAR contains the SUSY and soft-SUSY-breaking parameters. An extention 
of the SLHA conventions is needed here. The new entries are: 



61 for A 

62 for k 

63 for A\ 

64 for A K 

65 for yU c fj = A (S) 



Note that in NMHDECAY 1.0, the switch 23 (corresponding to the MSSM fi parameter) 
was used for n c g. Introducing a new switch for this parameter allows, in principle, non 
zero values for both fi and /i e ff , although such a scenario is not implemented in NMHDE- 
CAY (the main motivation for the NMSSM is to get rid of the MSSM // parameter). A 
non zero value for ji in switch 23 is simply ignored. 

The output files of nmhdecay_slha are spectr.dat, decay.dat and omega.dat. The direc- 
tory NMHDECAY contains test versions of these files - spectr.dat. test (see Table 2 below), 
decay.dat. test and omega.dat. test - corresponding to the sample input file slhainp.dat. The 
content of spectr.dat is as follows: 

a) BLOCK SPINFO is followed by warnings (switch 3) if any phenomenological constraint 
is violated. This segment of the output also displays error messages (switch 4) if any of 
the Higgs, squark or slepton states have a negative mass squared. No spectrum output is 
produced in this latter case. 

b) BLOCK SMINPUTS contains a printout of the Standard Model input parameters. The 
numerical values for Mw, m s , m c , V us , V c b and V u b, that have no SLHA numbers, also 
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appear in lines subsequent to # SMINPUTS Beyond SLHA. 

c) BLOCK M IN PAR is followed by a printout of the value of tan/3. 

d) BLOCK EXTPAR displays a printout of the SUSY and soft-SUSY-breaking parameters. 

e) BLOCK MASS contains the masses of all Higgs and sparticle states. There, one finds 
several essential NMSSM generalizations of the SLHA conventions. The new entries, with 
proposed PDG codes, are 

45 for the third CP-even Higgs boson, 

46 for the second CP-odd Higgs boson, 
1000045 for the fifth neutralino. 

f) BLOCK LOWEN contains observables relevant for precision experiments at low energy, 
at present only the computed value for BR(6 — > 57) (switch 1). 

g) The Higgs mixings in the CP-even sector follow BLOCK NMHMIX and those in the CP- 
odd sector follow BLOCK NMAMIX. Both segments are required in order to parameterize 
the mixing in the enlarged Higgs sector in the NMSSM. The meaning of the matrix 
elements (i, j = 1, 2, 3) and (i, j = 1, 2, 3) is as follows. 

• According to the SLHA conventions the Higgs weak eigenstates H u , H d are denoted 
by H 2 , Hi, respectively. (Inside the Fortran code the Higgs states H u , H d are 
denoted by Hi, H 2 , which is of no relevance for the SLHA output.) Hence, for the 
purpose of the SLHA output, the CP-even Higgs states are numbered by Sf eak = 
(H dR , H uR , S R ) (R refers to the real component of the field). If hi are the mass 
eigenstates (ordered in mass), the convention is hi = SijSJ eak . 

• In the CP-odd sector the weak eigenstates are H u j, H d i, Si (I for imaginary com- 
ponent). Again, for the purpose of the SLHA output (NEW according to SLHA2 
[21]!), the CP-odd Higgs states are denoted by H uI = H 2 i, H dI = Hu and Si = H31. 
The mass eigenstates are ai where a± is the Goldstone mode G, and the two phys- 
ical states a 2 , (13 are ordered in mass. Then the elements of P^- are defined as 

(X{ PijH ji . 
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h) BLOCK STOPMIX, SBOTMIX and STAUMIX contain the mixing matrices of the stop 
squarks, sbottom squarks and stau sleptons respectively, as defined in the MSSM. 

i) BLOCK NMNMIX is followed by a printout of the obvious generalization of the 4x4 
MSSM neutralino mixing matrix to the 5x5 NMSSM neutralino mixing matrix (with 
real entries); BLOCK UMIX and BLOCK VMIX are followed by printouts of the U and V 
matrices as defined in the MSSM. 

h) As indicated, the subsequent BLOCKs GAUGE, YU, YD, YE, L/K, AU, AD, AE, AL/AK 

and MSOFT contain the couplings and soft SUSY breaking terms at the SUSY breaking 
scale, and the following BLOCKs with suffix GUT the corresponding parameters at the 
GUT scale. 

The output file decay.dat gives the decay widths of all Higgs states into two particles, 
using the SLHA conventions and the above generalizations of the PDG codes both for the 
decaying particle and the final state particles. BLOCK DCINFO gives informations about 
the decay package (NMHDECAY version 2.1). 

The output file omega.dat contains the relic density of the lightest neutralino (if 
0MGFLAG=1). BLOCK RDINFO gives informations about the relic density package (Mi- 
crOMEGAs version 1.3). BLOCK RELDEN is then followed by switch 1 corresponding to 
the relic density Qh 2 and switch 2 corresponding to a warning in case the relic density of 
dark matter could not be computed or is excluded by WMAP bounds. 

3.2 NMHDECAY_SCAN 

The program n m h decay _scan uses the input file scaninp.dat, a version of which (scan- 
inp.dat.test) is downloaded automatically with the Fortran code (see Table 3 below). In 
this input file, the following parameters must be specified: 

• the total number of points to be scanned in parameter space; 

• the output format flag is for "short", corresponding to simple rows of numbers 
per allowed point in parameter space, and 1 for "long" , as described below; 
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• the flag OMGFLAG (0 for no relic density computation, 1 for relic density computa- 
tion using MicrOMEGAs); 

• lower and upper limits for the NMSSM parameters A, k, tan/3, fj, e s, A x and A K ; 

• the soft squark and slepton masses, trilinear couplings and the gaugino masses over 
all of which no scan is performed. 

The scan in parameter space uses a random number generator, such that all NMSSM 
parameters are randomly chosen point by point in the parameter space within the spec- 
ified limits. The standard model parameters (a s (Mz), Gf, ck^, the lepton masses m T 
and m^, the gauge boson masses Mz and Mw, the quark pole masses m s , m c , and m t , 
the running bottom quark mass mbirrib) and the CKM matrix elements V us , V c b and V u b) 
are specified in the subroutine INITIALIZE. 

The output file containing the physical parameters is always called scanout.dat, re- 
gardless of the output format chosen. The numbers printed out for the output format 
(recommended for scans over more than 10 points in parameter space) should be edited 
according to the user's needs (see the section of the program in the subroutine OUTPUT 
following the comment line The following output can be edited according to the user's needs). 
The output format 1 is easily readable and shows 

• the NMSSM parameters for each point as used as input, at the SUSY breaking scale, 
and at the GUT scale; 

• possible warnings in case any phenomenological constraint is violated, or error mes- 
sages ("fatal" errors) in case any of the Higgs or sfermion states has a negative mass 
squared (in which case no additional output is produced); 

• for each of the six Higgs states, their mass, their decomposition in the basis of 
interaction eigenstates (H u ,Hd,S), their reduced couplings to gauge bosons (CV), 
up type quarks (CU), down type quarks (CD), two gluons (CG) and two photons 
(CGA) (all relative to a standard model Higgs boson with the same mass), their 
branching ratios (where "Higgses" denote all possible two Higgs final states, and 
"sparticles" all possible two particle neutralino/chargino/sfermion final states), and 
their total width; 
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• the neutralino, chargino and gluino masses (all masses in GeV), as well as the 
neutralino composition in the basis ip Q = (—iX 1 , — i\ 2 , ip^, ip s )] 

• the pole masses and cosine of the mixing angle for the stop and sbottom states, as 
well as the pole masses of the squarks of the first two generations; 

• the masses and cosine of the mixing angle for the stau states, as well as the mass of 
the tau sneutrino and the masses of the first two generation sleptons; 

• the dark matter relic density (if 0MGFLAG=1); 

• the branching ratio BR(6 — > S7). 

The output file scanerr.dat shows how many of the points in parameter space have 
avoided fatal errors or violations of phenomenological constraints, and the range in the 
NMSSM parameter space over which points have passed all these tests. 

Users who wish to call a subroutine as a function of the Higgs or sparticle properties 
(masses, mixing angles and other quantities computed during the course of the scan) 
should use the parameters and common blocks found in the subroutine OUTPUT. The 
comments at the beginning of the main program should allow easy identification of all 
the parameters, branching ratios, mixing angles and so forth that would be of potential 
interest for inputting into a user's subroutine. 

4 Summary and outlook 

The versions 2.0+ of NMHDECAY are both a calculator of the NMSSM Higgs and spar- 
ticle spectrum, and a Higgs decay package. A quite unique feature of NMHDECAY is the 
check of each point in parameter space against limits on Higgs bosons and sparticles from 
accelerator experiments (LEP and Tevatron), that would have been impossible without 
the help of numerous collegues participating in these experiments. 

Apart from additional radiative corrections to the Higgs masses, the new features of the 
version 2.0 of NMHDECAY are the inclusion of Higgs decays into squarks and sleptons, 
the computation of the squark, gluino and slepton masses and mixings, the computation 
of the dark matter relic density and the branching ratio BR(6 — > S7). 
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Whereas the computation of the dark matter relic density is quite reliable (as in the 
version 1.3 of MicrOMEGAs), the result for the branching ratio BR(6 — > 57) has still to 
be interpreted with care, notably for larger (> 5) values of tan (3. 

It is clear that further improvements of NMHDECAY would be desirable: more higher 
order corrections to both the Higgs masses and decay widths would be welcome, with the 
aim to reach the present accuracy in the MSSM. (Tests of NMHDECAY in the MSSM 
limit A, k — > with /i e s fixed indicate, however, that the deviation of the mass of the 
lightest CP even Higgs boson w.r.t. corresponding MSSM calculations, for the same CP 
odd Higgs pole mass and sparticle spectra, is limited to about 3% and mostly much 
smaller.) 

Informations on low energy precision observables are included only in the form of a 
rough calculation of BR(b — > S7) , which should certainly be improved. Also the anomalous 
magnetic moment of the muon as well as Ap should be computed. 

We plan to treat these issues in the near future. 

Finally, it would be useful to be able to choose (universal) soft SUSY breaking terms 
at a GUT scale as in a mSUGRA version of the NMSSM. A corresponding code NMSPEC 
is in preparation. 
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# INPUT FILE FOR NMHDECAY VERSION 2.1 

# BASED ON SUSY LES HOUCHES ACCORD II 
BLOCK MODSEL 

3 1 # NMSSM PARTICLE CONTENT 

9 # FLAG FOR CALL OF MICROMEGAS (0=N0, 1=YES) 

BLOCK SMINPUTS 



1 


1. 


27920000E+02 


# 


ALPHA_EM~-1(MZ) 


2 


1, 


. 16639000E-05 


# 


GF 


3 


1. 


. 17200000E-01 


# 


ALPHA_S(MZ) 


4 


9. 


. 11870000E+01 


# 


MZ 


5 


4. 


21400000E+00 


# 


MB (MB) 


6 


1, 


75000000E+02 


# 


MTOP (POLE MASS) 


7 


1. 


77710000E+00 


# 


MTAU 


: MINPAR 











4. 


68008177E+02 


# 


REN. SCALE 


3 


3. 


OOOOOOOOE+00 


# 


TANBETA 



BLOCK EXTPAR 



1 


8, 


. 16527293E+01 


# 


Ml 


2 


1. 


53892994E+02 


# 


M2 


3 


4. 


72879088E+02 


# 


M3 


11 


-3. 


96388267E+02 


# 


ATOP 


12 


-7, 


44522774E+02 


# 


ABOT 


13 


-3. 


. 22488355E+02 


# 


ATAU 


65 


3, 


.53337956E+02 


# 


MU AT MZ 


31 


2. 


40063601E+02 


# 


M_eL 


32 


2, 


40063601E+02 


# 


M_muL 


33 


2. 


, 39934370E+02 


# 


M_tauL 


34 


2 


.13138181E+02 


# 


M_eR 


35 


2. 


, 13138181E+02 


# 


M_muR 


36 


2. 


. 12838756E+02 


# 


M_tauR 


41 


4. 


.74810720E+02 


# 


M_qlL 


42 


4. 


.74810720E+02 


# 


M_q2L 


43 


4. 


.20057207E+02 


# 


M_q3L 


44 


4. 


61771920E+02 


# 


M_uR 


45 


4. 


.61771920E+02 


# 


M_cR 


46 


3, 


36293030E+02 


# 


M_tR 


47 


4. 


60437692E+02 


# 


M_dR 


48 


4. 


.60437692E+02 


# 


M_sR 


49 


4. 


60051305E+02 


# 


M_bR 


61 


2. 


.OOOOOOOOE-01 


# 


LAMBDA AT MZ 


62 


1. 


.47434810E-01 


# 


KAPPA AT MZ 


63 


-7. 


50350192E+01 


# 


A_LAMBDA AT MZ 


64 


-2. 


. 18338774E+00 


# 


A_KAPPA AT MZ 



Table 1: Sample slhainp.dat file. 
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# NMHDECAY OUTPUT IN SLHA FORMAT 

# Info about spectrum calculator 



BLOCK SPINFO 

1 NMHDECAY 

2 2.1 

# Input parameters 
BLOCK MODSEL 

3 1 
BLOCK SMINPUTS 



# Program information 

# spectrum calculator 

# version number 



# NMSSM PARTICLE CONTENT 



1 


1 


27920000E+02 


# 


ALPHA_EM~-1(MZ) 


2 


1 


16639000E-05 


# 


GF 


3 


1 


17200000E-01 


# 


ALPHA_S(MZ) 


4 


9 


11870000E+01 


# 


MZ 


5 


4 


21400000E+00 


# 


MB (MB) 


6 


1 


75000000E+02 


# 


MTOP (POLE MASS) 


7 


1 


77710000E+00 


# 


MTAU 



# SMINPUTS Beyond SLHA: 

# MW: 0.80420000E+02 

# MS: 0.19000000E+00 

# MC: 0.14000000E+01 

# VUS: 0.22000000E+00 

# VCB: 0.40000000E-01 

# VUB: 0.40000000E-02 
BLOCK MINPAR 

3 3.00000000E+00 
BLOCK EXTPAR 

1 8.16527293E+01 

2 1.53892994E+02 

3 4.72879088E+02 

11 -3.96388267E+02 

12 -7.44522774E+02 

13 -3.22488355E+02 
65 3.53337956E+02 

31 2.40063601E+02 

32 2.40063601E+02 

33 2.39934370E+02 

34 2.13138181E+02 

35 2.13138181E+02 

36 2.12838756E+02 

41 4.74810720E+02 

42 4.74810720E+02 

43 4.20057207E+02 

44 4.61771920E+02 

45 4.61771920E+02 

46 3.36293030E+02 

47 4.60437692E+02 

48 4.60437692E+02 

49 4.60051305E+02 

61 2.00000000E-01 

62 1.47434810E-01 

63 -7.50350192E+01 

64 -2.18338774E+00 



# TANBETA 



Ml 
M2 
M3 
ATOP 

# ABOTTOM 

# ATAU 

# MU 

# LEFT SELECTRON 

# LEFT SMUON 

# LEFT STAU 

# RIGHT SELECTRON 
RIGHT SMUON 
RIGHT STAU 
LEFT 1ST GEN. 
LEFT 2ND GEN. 
LEFT 3RD GEN. 
RIGHT U-SQUARKS 
RIGHT C-SQUARKS 

# RIGHT T-SqUARKS 

# RIGHT D-SqUARKS 
RIGHT S-SqUARKS 
RIGHT B-SqUARKS 
LAMBDA 
KAPPA 
A_LAMBDA 
A KAPPA 



SqUARKS 
SqUARKS 
SqUARKS 
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BLOCK MASS # Mass spectrum 
# PDG Code mass 



25 


9. 


74781199E+01 


35 


4. 


48583968E+02 


45 


5. 


. 16024735E+02 


36 


2. 


19434909E+01 


46 


4. 


50932452E+02 


37 


4. 


55175450E+02 


1000001 


4. 


96428647E+02 


2000001 


4. 


80320535E+02 


1000002 


4, 


91385003E+02 


2000002 


4. 


80115856E+02 


1000003 


4, 


. 96428647E+02 


2000003 


4, 


80320535E+02 


1000004 


4. 


91385003E+02 


2000004 


4, 


80115856E+02 


1000005 


4, 


40278 168E+02 


2000005 


4, 


. 77873225E+02 


1000006 


3, 


. 03200260E+02 


2000006 


5. 


, 24782 130E+02 


1000011 


2. 


. 437475 16E+02 


2000011 


2. 


16456923E+02 


1000012 


2. 


33286301E+02 


1000013 


2, 


, 437475 16E+02 


2000013 


2. 


. 16456923E+02 


1000014 


2, 


33286301E+02 


1000015 


2, 


. 15119676E+02 


2000015 


2 


. 445541 12E+02 


1000016 


2. 


33153313E+02 


1000021 


5. 


01506590E+02 


1000022 


7. 


55014082E+01 


1000023 


1. 


38801726E+02 


1000025 


-3. 


.58593816E+02 


1000035 


3. 


. 80957538E+02 


1000045 


5. 


. 24690943E+02 


1000024 


1. 


. 37945343E+02 


1000037 


3, 


79628706E+02 



# Low energy observables 
BLOCK L0WEN 

1 3.33858054E+00 

# 3*3 Higgs mixing 
BLOCK NMHMIX 



1 


1 


-3 


.33228276E- 


•01 


# 


1 


2 


-9 


.41055786E- 


■01 


# 


1 


3 


5. 


,80768815E- 


02 


# 


2 


1 


-9. 


. 29046853E- 


■01 


# 


2 


2 


3. 


, 17227652E- 


•01 


# 


2 


3 


-1. 


. 90364284E- 


■01 


# 


3 


1 


-1. 


.60719818E- 


■01 


# 


3 


2 


1. 


. 17390906E- 


•01 


# 


3 


3 


9. 


.79994140E- 


■01 


# 



particle 

# lightest neutral scalar 

# second neutral scalar 

# third neutral scalar 

# lightest pseudoscalar 

# second pseudoscalar 

# charged Higgs 

# ~d_L 

# ~d_R 

# ~u_L 

# ~u_R 

# ~s_L 

# ~s_R 

# ~c_L 

# ~c_R 

# ~b_l 

# ~b_2 

# ~t_l 

# ~t_2 

# ~e_L 

# ~e_R 

# ~nue_L 

# ~mu_L 

# ~mu_R 

# ~numu_L 

# ~tau_l 

# ~tau_2 

# ~nutau_L 

# ~g 

# neutralino (1) 

# neutralino (2) 

# neutralino (3) 

# neutralino (4) 

# neutralino (5) 

# chargino(l) 

# chargino(2) 



# BR(b -> s gamma) *10~4 



S_(l,l) 
S_(l,2) 
S_(l,3) 
S_(2,l) 
S_(2,2) 
S_(2,3) 
S_(3,l) 
S_(3,2) 
S_(3,3) 
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# 3*3 Pseudoscalar Higgs mixing 
BLOCK NMAMIX 



P_(l,l) 
P_(l,2) 
P_(l,3) 
P_(2,l) 

# P_(2,2) 

# P_(2,3) 

# P_(3,l) 

# P_(3,2) 

# P_(3,3) 



1 1 3.16227766E-01 
1 2 -9.48683298E-01 

1 3 0.00000000E+00 

2 1 -9.84143132E-02 
2 2 -3.28047711E-02 

2 3 -9.94604680E-01 

3 1 9.43564848E-01 
3 2 3.14521616E-01 
3 3 -1.03737795E-01 

# 3rd generation sfermion mixing 
BLOCK STOPMIX # Stop mixing matrix 

1 1 5.65184261E-01 # Rst_(l,l) 

1 2 8.24964697E-01 # Rst_(l,2) 

2 1 -8.24964697E-01 # Rst_(2,l) 
2 2 5.65184261E-01 # Rst_(2,2) 

BLOCK SBOTMIX # Sbottom mixing matrix 
1 1 9.94871315E-01 # Rsb_(l,l) 

1 2 1.01148731E-01 # Rsb_(l,2) 

2 1 -1.01148731E-01 # Rsb_(2,l) 
2 2 9.94871315E-01 # Rsb_(2,2) 

BLOCK STAUMIX # Stau mixing matrix 
1 1 1.82923365E-01 # Rsl_(l,l) 

1 2 9.83127175E-01 # Rsl_(l,2) 

2 1 -9.83127175E-01 # Rsl_(2,l) 
2 2 1.82923365E-01 # Rsl_(2,2) 

# Gaugino-Higgsino mixing 

BLOCK NMNMIX # 5*5 Neutralino Mixing Matrix 



1 


1 


9 


76939885E- 


01 


# 


N_ 


(1 


1 


1 


2 


-1 


21822079E- 


01 


# 


N_ 


(1 


2 


1 


3 


1 


55584479E- 


01 


# 


N_ 


(1 


3 


1 


4 


-8 


03150977E- 


02 


# 


N_ 


(1 


4 


1 


5 


9 


52878280E- 


03 


# 


N_ 


(1 


5 


2 


1 


1 


69572219E- 


01 


# 


N_ 


(2 


1 


2 


2 


9 


41177948E- 


01 


# 


N_ 


(2 


2 


2 


3 


-2 


44926790E- 


01 


# 


N_ 


(2 


3 


2 


4 


1 


58638593E- 


01 


# 


N_ 


(2 


4 


2 


5 


-1 


65528778E- 


•02 


# 


M_ 


(2 


5 


3 


1 


-4 


43118962E- 


02 


# 


N_ 


(3 


1 


3 


2 


7 


05108082E- 


02 


# 


M_ 


(3 


2 


3 


3 


6 


98006382E- 


01 


# 


N_ 


(3 


3 


3 


4 


7 


10367686E- 


01 


# 


N_ 


(3 


4 


3 


5 


3 


50645553E- 


02 


# 


N_ 


(3 


5 


4 


1 


-1 


21736484E- 


01 


# 


N_ 


(4 


1 


4 


2 


3 


06778393E- 


01 


# 


N_ 


(4 


2 


4 


3 


6 


47820106E- 


01 


# 


M_ 


(4 


3 


4 


4 


-6 


79460614E- 


01 


# 


N_ 


(4 


4 


4 


5 


9 


86388891E- 


02 


# 


N_ 


(4 


5 


5 


1 


7 


09990264E- 


03 


# 


N_ 


(5 


1 


5 


2 


-1 


60840169E- 


02 


# 


N_ 


(5 


2 


5 


3 


-9 


444860 10E- 


02 


# 


M_ 


(5 


3 


5 


4 


4 


57635935E- 


02 


# 


N_ 


(5 


4 


5 


5 


9 


94321905E- 


01 


# 


N_ 


(5 


5 
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BLOCK UMIX # Chargino U Mixing Matrix 



1 


1 


9 


.25480029E- 


-01 


# 


U_ 


(1,1) 


1 


2 


-3 


. 78796404E- 


-01 


# 


u_ 


(1,2) 


2 


1 


3 


.78796404E- 


-01 


# 


u 


(2,1) 


2 


2 


9 


.25480029E- 


-01 


# 


u 


(2,2) 


,0CK VMIX 


# Chargino 


V Mixing 


Matrix 


1 


1 


9 


.69142290E- 


-01 


# 


v_ 


(1,1) 


1 


2 


-2 


.46501971E- 


-01 


# 


v_ 


(1,2) 


2 


1 


2 


.46501971E- 


-01 


# 


v_ 


(2,1) 


2 


2 


9 


.69142290E- 


-01 


# 


v_ 


(2,2) 



# 



# GAUGE AND YUKAWA COUPLINGS AT THE SUSY SCALE 
BLOCK GAUGE 0= 4 . 68008177E+02 # (SUSY SCALE) 

1 3.60690129E-01 # gl(Q,DR_bar) 

2 6.45231920E-01 # g2(Q,DR_bar) 

3 1.10281535E+00 # g3(Q,DR_bar) 
BLOCK YU 0= 4.68008177E+02 # (SUSY SCALE) 

3 3 9.40962309E-01 # HTOP(Q,DR_bar) 
BLOCK YD 0= 4.68008177E+02 # (SUSY SCALE) 

3 3 4.53102195E-02 # HBOT(Q,DR_bar) 
BLOCK YE 0= 4.68008177E+02 # (SUSY SCALE) 

3 3 3.17949988E-02 # HTAU(Q,DR_bar) 
BLOCK L/K Q= 4.68008177E+02 # (SUSY SCALE) 

1 2.00784168E-01 # LAMBDA (Q,DR_bar) 

2 1.48000459E-01 # KAPPA (Q,DR_bar) 

# 

# SOFT TRILINEAR COUPLINGS AT THE SUSY SCALE 
BLOCK AU 0= 4.68008177E+02 # (SUSY SCALE) 

3 3 -3.96388267E+02 # ATOP 
BLOCK AD 0= 4.68008177E+02 # (SUSY SCALE) 

3 3 -7.44522774E+02 # ABOT 
BLOCK AE 0= 4.68008177E+02 # (SUSY SCALE) 

3 3 -3.22488355E+02 # ATAU 
BLOCK AL/AK Q= 4 . 68008177E+02 # (SUSY SCALE) 



1 


-8. 


.42676828E+01 


# 


ALAMBDA 


2 


-2 


.56530300E+00 


# 


AKAPPA 


# 

# SOFT MASSES 


AT 


THE SUSY SCALE 






BLOCK MSOFT Q= 


= 4.68008177E+02 


# 


(SUSY SCALE) 


1 


8. 


16527293E+01 


# 


Ml 


2 


1. 


53892994E+02 


# 


M2 


3 


4. 


72879088E+02 


# 


M3 


21 


5, 


81854464E+04 


# 


M_HD~2 


22 


-1 


.04142570E+05 


# 


M_HU~2 


23 


-1. 


. 32549686E+05 


# 


M_S~2 


31 


2. 


40063601E+02 


# 


M_eL 


32 


2. 


40063601E+02 


# 


M_muL 


33 


2. 


39934370E+02 


# 


M_tauL 


34 


2. 


. 13138181E+02 


# 


M_eR 


35 


2. 


13138181E+02 


# 


M_muR 


36 


2, 


12838756E+02 


# 


M_tauR 


41 


4. 


74810720E+02 


# 


M_qlL 


42 


4. 


74810720E+02 


# 


M_q2L 


43 


4. 


20057207E+02 


# 


M_q3L 


44 


4. 


61771920E+02 


# 


M_uR 


45 


4. 


61771920E+02 


# 


M_cR 


46 


3. 


36293030E+02 


# 


M_tR 


47 


4. 


60437692E+02 


# 


M_dR 


48 


4. 


60437692E+02 


# 


M sR 


49 


4. 


60051305E+02 


# 


M_bR 
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# MU_EFF: 3 . 53749043E+02 # (AT THE SUSY SCALE) 
# 

# GAUGE AND YUKAWA COUPLINGS AT THE GUT SCALE 
BLOCK GAUGEGUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 

1 7. 18934347E-01 # gl (MGUT,DR_bar) , GUT normalization 

2 7. 18934334E-01 # g2 (MGUT , DR_bar) 

3 7.09148313E-01 # g3 (MGUT , DR_bar) 
BLOCK YUGUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 

3 3 6.29869150E-01 # HTOP (MGUT , DR_bar) 
BLOCK YD GUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 

3 3 1.69658275E-02 # HBOT (MGUT , DR_bar) 

BLOCK YEGUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 

3 3 2. 17473793E-02 # HTAU (MGUT , DR_bar ) 
BLOCK L GUT /K GUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 

1 2.22229350E-01 # LAMBDA (MGUT, DR_bar) 

2 1.61430186E-01 # KAPPA (MGUT, DR_bar) 

# 

# SOFT TRILINEAR COUPLINGS AT THE GUT SCALE 
BLOCK AUGUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 

3 3 -2.00014623E+02 # ATOP 
BLOCK AD GUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 

3 3 -1.99980396E+02 # ABOT 
BLOCK AEGUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 

3 3 -1.99994867E+02 # ATAU 
BLOCK ALGUT/AKGUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 

1 -1.99976861E+02 # ALAMBDA 

2 -2.29490634E+01 # AKAPPA 

# 

# SOFT MASSES SQUARED AT THE GUT SCALE 

BLOCK MSOFTGUT MGUT= 2 . 39561249E+16 # (GUT SCALE) 



1 


2. 


.00004409E+02 


# 


Ml 


2 


2 


.00001910E+02 


# 


M2 


3 


1. 


.99989075E+02 


# 


M3 


21 


3. 


.99907796E+04 


# 


M_HD 


22 


4. 


.00122196E+04 


# 


M_HU 


23 


-1. 


.42907310E+05 


# 


M_S 


31 


3. 


.99987978E+04 


# 


M_eL 


32 


3. 


.99987978E+04 


# 


M_muL 


33 


3. 


.99987958E+04 


# 


M_tauL 


34 


4. 


.00004265E+04 


# 


M_eR 


35 


4. 


.00004265E+04 


# 


M_muR 


36 


4. 


.00004222E+04 


# 


M_tauR 


41 


3. 


.99950362E+04 


# 


M_qlL 


42 


3. 


.99950362E+04 


# 


M q2L 


43 


3. 


.99997864E+04 


# 


M_q3L 


44 


3. 


.99948593E+04 


# 


M_uR 


45 


3. 


.99948593E+04 


# 


M_cR 


46 


4. 


.00046814E+04 


# 


M_tR 


47 


3. 


.99956655E+04 


# 


M_dR 


48 


3. 


.99956655E+04 


# 


M_sR 


49 


3. 


.99956856E+04 


# 


M_bR 



Table 2: Corresponding spectr.dat. test output file. 
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# 

# Total number of points scanned 
# 

1 
# 

# Output format 0=short l=long (not recommended for big scannings) 
# 

1 
# 

# Computation of relic density using MicrOmegas (0=no, l=yes) 
# 

1 
# 

# lambda 
# 

.2D0 
.2D0 
# 

# kappa 
# 

. 147D0 
. 147D0 
# 

# tan(beta) 
# 

3. DO 
3. DO 
# 

# mu 
# 

353. DO 
353. DO 
# 

# A_lambda 
# 

-75. DO 
-75. DO 
# 

# A_kappa 
# 

-2.2D0 
-2.2D0 
# 

# Remaining soft terms (no scan) 
# 

mQ3= 420. DO 
mU3= 336 . 3D0 
mD3= 460. DO 
mL3= 240. DO 
mE3= 213. DO 
AU3= -396. 4D0 
AD3= -744. 5D0 
AE3= -322. 5D0 
mQ= 474. 8D0 
mU= 461.8D0 
mD= 460. 4D0 
mL= 240. DO 
mE= 213. 1D0 
Ml= 81.7D0 
M2= 153. 9D0 
M3= 472. 9D0 



Table 3: The scaninp.dat file for sample parameter scan. 
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